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ABSTRACT 

As a neutron star spins down, the nuclear matter continuously is converted into quark 
matter due to the core density increase and then latent heat is released. We have inves- 
tigated the thermal evolution of neutron stars undergoing such deconfinement phase 
transition. We have taken into account the conversion in the frame of the general 
theory of relativity. The released energy has been estimated as a function of change 
rate of deconfinement baryon number. Numerical solutions to cooling equation are 
obtained to be very different from the without heating effect. The results show that 
neutron stars may be heated to higher temperature which is well-matched with pul- 
sar's data despite onset of fast cooling in neutron stars with quark matter core. It is 
also found that heating effect has magnetic field strength dependence. This feature 
could be particularly interesting for high temperature of low-field millisecond pul- 
sar at late stage. The high temperature could fit the observed temperature for PSR 
J0437-4715. 
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1 INTRODUCTION 

The composition of Neutron star (NS) interior is still poorly known due to the uncertainties of nu- 
clear physics. NS cooling is an important tool for the study of dense matter. By comparing cooling 
models with thermal emission data from observations, we can gain insight into the equation of state 
(EOS) of dense matter inside NSs. After a supernova explosion a newly formed NS first cools via 
various neutrino emission mechanisms before the surface photon radiation take over. Only slower 
neutrino cooling mechanisms such as the modified Urea, plasmon neutrino and bremsstrahlung 



* E-mail:kangmiao@phy.ccnu. edu.cn 
t E-mail:zhxp@phy.ccnu. edu.cn 



2 Miao Kang , Xiao-Ping Zheng 

processes, should occur when the interior density is not high. However, for higher core density the 
more 'exotic' extremely fast cooling processes take over. Cooling of hybrid stars (NSs with quark 
matter core) is just one of those examples. 

However, heating effect on the cooling of compact stars is an important factor. Several heating 
mechanisms have been extensively discussed, for example, the dissipation of rotational energy due 
to viscous damping (Zheng & Yu ( 2006)) and mutual friction between superfluid and normal com- 
ponents of the star(Shibazaki & Lamb (1989)) and release of strain energy stored by the solid crust 
due to spin-down deformation(Cheng et al. (|19921) ). These heating processes are closely related 
with the rotation evolution of a star. It is well-known that a neutron star spin-down due to magnetic 
dipole radiation. Because of spin-down compression, the interior density of the star will gradually 
increase. For hybrid stars, this results in, little by little, the transformation of hadron matter into 
quark matter in the interior. It will lead to the release of latent heat if the transition is the first-order 
one. The generation of the energy increases internal energy of the star. It will be called decon- 
finement heating (DH). Deconfinement process has ever been investigated in strange stars, where 
neutron drops at bottom of a crust drip on quark matter surface to be instantaneously dissolved 
into quark matter (Haensel & Zdunik (fT99~Ti Yuan & Zhang (fT999l) . Yu & Zheng (120061) ). What 
we will here deliberate is slow transition process: nuclear matter undergoing from hadron matter 
phase (HP) to mixed hadron-quark matter phase (MP) and then to quark matter phase (QP) with 
gradually increasing density. 

From the point of view above, the change of the internal structure of the compact star due 
to rotation has to be evaluated within general relativity theory. Using the method of perturba- 
tion theory(Hartle (119671) , Chubarian et al (2000)), we are going to investigate evolution of phase 
transition region during the spin-down evolution of the rotating star. And then we also calculate 
changes in confined baryon numbers and quark numbers. The heat luminosity can be estimated 
with proportion to the change rates. 

Of course, the star is not quite in weak interaction equilibrium state during spin down. The 
departure from the chemical equilibrium indeed lead to the rotochemical heating in a rotating 
neutron star (Reisenegger (1995), (2006)), as well as a hybrid star. The direct Urea processes 
would be triggered and dominates the rotochemical heating in our models. We find, however, the 
effect of rotochemical heating is much smaller than the deconfinement heating by estimating the 
both net heating rates for the sequences of different parameters. Hence, we don't consider the 
effect of rotochemical heating in following study. 

In this work, we take Glendenning's hybrid stars model (Glendenning (1997)). The hadronic 
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matter equations of state(EOS) in the framework of the relativistic mean field theory and MIT 
bag model of quark matter are used to construct the model of hybrid stars, but medium effect of 
quark matter has been considered in quasiparticle description(Schertle et al (119971) ). We choose 
the simplest possible nuclear matter composition, namely neutrons, protons, electrons, and muons 
(npe/i matter) and ignore superfluidity and superconductivity. 

The plan of this paper is as follows. In Sec. 2 we introduce rotating hybrid star model. The DH 
effect is considered in Sec. 3. The cooling curves and the corresponding explanations are presented 
in Sec. 4. The conclusion and discussions are summarized in Sec. 5. 



2. ROTATING HYBRID STARS MODEL 

Quark deconfiment phase transition is expected to occur in neutron matter at densities above the 
nuclear saturation density n b = 0.16/ra -3 . Since many theoretical calculations have suggested 
that deconfinement transition should be of first order in low-temperature and high-density area 
(Pisalski & Wilczek (119841) and Gavai et al. (119871) ). one may expect a MP during the transition. 
Most of the approaches to deconfinement matter in neutron star matter use a standard two-phase 
description of EOS where the HP and the QP are modelled separately and resulting EOS of the 
MP is obtained by imposing Gibbs conditions for phase equilibrium with the constraint that baryon 
number as well as electric charge of the system are conserved (Glendenning (119921) . (119971) ). 

The Gibbs condition for mechanical and chemical equilibrium at zero temperature between the 
HP and the QP reads 

PHp(Mn,Ve) = P QP (Mn,Ve)- (1) 

where p HP is pressure of HP and p QP is the pressure of QP. We use the EOS of the relativistic 
mean field model (Glendenning (119971) ) for hadron matter and employ a effective mass bag-model 
EOS for quark matter(Schertle et al. (11997b ). Only two independent chemical potentials remain ac- 
cording to the corresponding two conserved charges of the /3-equilibrium system. The total baryon 
number N B as well as electrical charge Q 
N B 

n B = — = X"qp + (1 ~x)n H p- (2) 
Q 

= y =XQqp + (1 ~X)Qhp- (3) 

where x - Vq/V is the volume fraction of quark matter in the MP. Taking the charge neutral EOS 
of the HP, Eq.(l), (2) and (3) for MP and the charge neutral EOS of the QP, we can construct the 
full hybrid star EOS. In Fig.l we show the model EOS with deconfinement transition which is 
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the typical scheme of a first order transition at finite density with MP. We choose the parameters 
for hadronic matter EOS which have given by Glendenning (1997 ) and quark matter EOS with s 
quark mass m s = 150MeV, bag constant 5 1/4 = 160MeV, coupling constant g = 3. 

With the evaluated hybrid star EOS presented above we now turn to analyse the structure of the 
corresponding rotating hybrid stars. Using the Hartle's perturbation theory (119671) , Chubarian et al 
( 2000} have studied the change of the internal structure of the hybrid stars due to rotation. In this 
paper, we also apply Hartle's approach to investigate the structure of rotating hybrid stars. Hartle's 
formalism is based on treating a rotating star as a perturbation on a non-rotating star, expanding 
the metric of an axially symmetric rotating star in even powers of the angular velocity Q.. The 
metric of a slowly rotating star to second order in the angular velocity Q, can be written as 

ds 2 = -e^\l + 2(h + h 2 P 2 )]dt 2 + e*»[l + ^ + ™f* V 2 

(r - 2M(r)) 

+r 2 [l + 2(v 2 - h 2 )P 2 ]{d6 2 + sin 2 9[d<f> - w(r, 8)dt] 2 } + 0(Q 3 ) (4) 

Here e v{r) ,e Mr) and M(r) are functions of r and describe the non-rotating star solution of the Tolman- 
Oppenheimer-Volkov (TOV) equations. P 2 = P 2 (9) is the 1 = 2 Legendre polynomials, oj is the 
angular velocity of the local inertial frame and is proportional to the star's angular velocity f2, 
whereas the perturbation functions h , h 2 , mo, m 2 , v 2 are proportional to Q 2 . we assume that matter 
in the star is described by a perfect fluid with energy momentum tensor 

T>" = (e + P)u M u v + Pg" v (5) 

The energy density and pressure of the fluid are affected by the rotation because the rotation 
deforms the star. In the interior of the star at given (r, 9), in a reference frame that is momentarily 
moving with the fluid, the pressure and energy density variation is respectively 

6P(r, 9) = [e(r) + P(r)] [p* + p* 2 P 2 (9)] (6) 

Se(r, 9) = ^[e(r) + P(r)][p*o + P*2 P 2(Q)] (7) 

here, pi and p* 2 are dimensionless functions of r, proportional to Q. 2 , which describe the pressure 
perturbation. The rotational perturbations of the star's structure are described by the functions 
ho, niQ,p* Q , h 2 , m 2 , v 2 ,p* 2 . These functions are calculated from Einstein's field equations. The effect 
of rotation described by the metric on the shape of the star can be divided into contributions: A 
spherical expansion which changes the radius of the star, and is described by the functions h 
and m . The other part is a quadrupole deformation, described by functions h 2 , v 2 and m 2 . As a 
consequence of these contributions, the difference between the gravitational mass of the rotating 
star and the non-rotating star with the same central pressure is 
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J 2 

6M grav = m (R) + — (8) 
The change in the radius of the star is given by 

8R = UR) + HR)Pi(.0) (9) 

We wish to study sequences of the rotating stars with constant total baryon number at variable spin 
frequency v = Q./2n. The expansion of total baryon numbers in powers of f2 is 

N B = N° + SN B + 0(Q 4 ) (10) 

where 

r R 

N° B = n B (r)[l - 2M(r)/r]- 1/2 4nr 2 dr (11) 
Jo 

is the total baryons number of non-rotating star and 



= -Lf(i- 

m N Jo 



2M(r) N _ 1/2fn m () (r) 1 _ 2 



8N B = — (1 — r 1/z {[l + ^- + -^[Q-o>(r)]V v ]mW) 

r-2M(r) 3 

^(^^(^ (12) 

here m/v is the rest mass per baryon. To construct constant baryon number sequences, we first solve 
the TOV equations to find the non-rotating configuration for giving a central pressure P(r=0). And 
then, for an assigned value of the angular velocity Q. the equations of star structure are solved to 
order Q. 2 , imposing that the correction to the pressure p* Q (r = 0) being not equal to zero. The value 
of p*(r = 0) is then changed until the same baryon number as that non-rotating star is obtained. 

The results for the stability of rotating hybrid star configurations with possible deconfinement 
phase transition according to the EOS described above are shown in Fig. 2, where the total gravita- 
tion mass is given as functions of the equatorial radius and the central baryon number density for 
static stars as well as for stars rotating with the maximum rotation frequency y^. The dotted lines 
connect configuration with the same total baryon number and it becomes apparent that the rotating 
configurations are less compact than the static ones. In order to explore the increase in central den- 
sity due to spin down, we create sequences of hybrid star models. Model in a particular sequence 
have the same constant baryon number, increasing central density and decreasing angular veloc- 
ity. Fig. 3 displays the central density of rotating hybrid stars with different gravitational mass at 
zero spin, as a function of its rotational frequency. In the interior of these stars, the matter can be 
gradually converted from the relatively incompressible nuclear matter phase to more compressible 
quark matter phase. 
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3. DECONFINEMENT HEATING 

As the star spins down, the centrifugal force decreases continuously, increasing its internal density. 
At this occurrence, the nuclear matter continuously converts into quark matter in an exothermic 
reaction, i.e. n — > u+2d, p —> 2u+d, s quarks immediately appear after weak decay. Fig. 3 identifies 
the fact that quarks are accumulating in the interior of the star with decreasing rotation frequency 
v. In a particular sequence with constant baryon number, we can calculate the deconfinement 
baryon number N q of the star with the varying rotation frequency. The expression of N q is similar 
to that the total baryon number Nb, except that the baryon number density n B is displaced with 
the deconfinement baryon number density n q . Deconfinement baryon number of 1.4M S star, for 
example, is plotted in Fig.4. The analytic expression fits as 

N q = iVj(l - 0.7 16vj + 0.055v^ - 0.032^) (13) 

where N® « 0.22N Q is the baryon number of quarks for the static configuration and V3 = v/10 3 Hz. 
We can also derive similar expressions for sequences of hybrid star models. Left panel in Fig.l 
indicates release of latent heat because there is the reduction in enthalpy per baryon at transition. 
Assumed the average value of release energy per nucleon that is transforming into quarks q n , the 
total latent heat per time can be written as 
dN a 

H dec {t) = q n —±v(t) (14) 
dv 

with 

y=-i^uV sin 2 (15) 

is induced by magnetic dipole radiation (MDR), where / is the stellar moment of inertia, /u = ^BR 3 
is the magnetic dipole moment, and 9 is the inclination angle between magnetic and rotational axes, 
q n can be estimated as the order of 0.1 MeV through contrast of HP enthalpy to corresponding MP 
enthalpy for same region of baryon number density, which is uniquely defined for given stellar 
mass in our work. 



4. COOLING CURVES 

The cooling is realized via two channels - by neutrino emission from the entire star body and 
by transport of heat from the internal layers to the surface resulting in the thermal emission of 
photons. Neutrino emission is generated in numerous reactions in the interiors of neutron stars, 
as reviewed, by Page et al. (2005 ). For the calculation of cooling of the hadron part of the hybrid 
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star we use the main processes which are nucleon direct Urea (NDU) and nucleon modified Urea 
(NMU) and nucleon bremsstrahlung (NB). If the proton and electron Fermi momenta are too 
small compared with the neutron Fermi momenta, the NDU process is forbidden because it is 
impossible to satisfy conservation of momentum. Under typical conditions one finds that the ratio 
of the number density of protons to that of nucleons must exceed about 0.1 1 for the process to be 
allowed. Medium effects and interactions among the particles modify this value only slightly but 
the presence of muons raise it to about 0.15. We apply the second condition in our calculations. 
For the calculation of cooling of the quark matter we consider the most efficient processes: the 
quark direct Urea (QDU) processes on unpaired quarks, the quark modified Urea (QMU) and the 
quark bremsstrahlung (QB). Nucleon superfluidity and quark superconductivity are not included 
in the model. Considering the energy equation of the star, the cooling equation can be written as 

C v — = -L v -L y +H (16) 
at 

where CV is the total stellar heat capacity, it incorporates neutron contribution, quark contribution 
and electron contribution. The term H indicates the heating energy per unit time, in our work 
H = H dec , L v is the neutrino luminosity, and L y is the surface photon luminosity given by 

L y = AnR 2 o-T A s , (17) 

here cr is the Stefan-Boltzmann constant and T s is the surface temperature. The surface temperature 
is related to internal temperature by a coefficient determined by the scattering processes occurring 
in the crust. We apply an formula which is demonstrated by Gudmundsson et al. (1983 ), 

T s = 3.08 x lO 6 ^^ 1 - 5495 (18) 

where g St u is the proper surface gravity of the star in units of I0 u cms~ 2 . The gravitational red-shift 
is also taken into account. Then the effective surface temperature detected by a distant observer is 
yoo _ y s j i _ R g /R, where R g is the gravitational stellar radius. In our calculation, we choose the 
initial temperature To = 10 9 k, and the magnetic tilt angle 6 = 45°. 

Fig. 5 shows the cooling behavior of a 1.4M Q hybrid star for different magnetic fields (10 9 - 
10 13 G). In this figure, q n is taken to be 0.1 MeV. It is evident that the DH increase the surface 
temperature dramatically. This is extremely different from fast cooling scenario (solid curve in 
Fig. 5). We find that there is a quite clear magnetic-field dependence of the curves. It is determined 
by the properties of MDR. Solutions to equation (15) show the magnetic-field dependence of 
the spin frequencies. The strong field strength induces a rapid spin-down at the beginning while 
the low field strength leads to only obvious spin-down at the older ages (see Fig. 3 in Zheng et 
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al (2006)). Hence, a large quantity of nuclear matter is deconfined into quark matter during the 
earlier and shorter period for strong-field case. On the contrary, the plateau can form at late period 
for weak-field case. 

In Fig. 6 we present the cooling behavior of different mass stars for B=10 12 with and without 
DH. The observational data, taken from tables 1 and 2 in Page et al. (l2004|) . have been shown in 
the figure. The theoretical curves are consistent with the observational data. The heating rate has a 
regular mass dependence. More massive neutron star have larger heating rate for the neutron stars 
gravitational mass M > 1.2M Q . 

Fig. 7 shows the evolution of the surface temperature of different mass stars for the weak 
(B=10 9 G) magnetic field. In the cases of weak field, stars could maintain high temperatures even 
at older ages (> 10 6 yrs). This feature may give an illustration of the inferred temperature for PSR 
J0437-4715 (Kargaltsev et al. (1200410 . although the quantitative analysis needs future study. 

5. CONCLUSIONS AND DISCUSSIONS 

The thermal evolution of rotating hybrid stars with DH have been investigated in this work. Using 
Hartle's perturbative approach, we have calculated the change of internal structure of rotating 
hybrid stars. The nuclear matter can continuously be converted into quark matter to release latent 
heat during the spins down of star. The heat luminosity can be estimated as proportion to the 
change rates of quark number. The results show the DH of rotating stars leads to good agreement 
with the observed data in the case of enhanced cooling (That is the onset of direct Urea process). 
We also found that heating effect has magnetic field strength dependence. For those stars with weak 
fields (< 10 10 G), our results show that they can maintain a high temperature (> 10 5 K) at older ages 
(t ~ I0 l0 yrs) such as PSR J0437-4715. It may be particularly interesting for high temperature of 
weak- field stars at late stage. 

We here calculate the released latent heat by regarding the heat release per nucleon as a pa- 
rameter. However, the heat release per nucleon should be a density-dependent quantity. Numerical 
calculations are necessary for heat luminosity. Precise fittings need consider improvement of the 
star model, the inclusions of such superfluidity, superconductivity and tension effect in MP. These 
will be the future works. 
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Figure 2. Gravitational mass M as a function of the equatorial radius (left figure) and the central density (right figure) for rotating hybrid stars 
configurations with a deconfinement phase phase transition. The solid curves correspond to static configurations, the dashed ones to those with 
maximum rotation frequency vjj. The lines between both extremal cases connect configurations with the same total baryon number. 
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Figure 4. The number of converting quark into baryon as a function of rotational frequency for 1 .4 M Q rotating hybrid star. 
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Figure 5. Cooling curves of 1.4 M Q hybrid star with DH for various magnetic fields (curve a: 10 13 , b: 10 12 , c: 10 11 , d: 10 10 , e: 10 9 ) and the curves 
without DH (solid curve). 




Figure 6. Cooling curves of neutron stars with DH for different star mass and B=10 12 G (upper curves) and the curves without DH (lower curves). 
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Figure 7. Cooling curves of neutron stars with DH for different star mass and B=10 5 G (upper curves) and the curves without DH (lower curves). 



